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Abstract. We discuss how well a given convex body B in a real d-dimensional vector 
space V can be approximated by a set X for which the membership question: "given 
an x G V, does x belong to X?" can be answered efficiently (in time polynomial 
in d). We discuss approximations of a convex body by an ellipsoid, by an algebraic 
hypersurface, by a projection of a polytope with a controlled number of facets, and by 
a section of the cone of positive semidefinite quadratic forms. We illustrate some of 
the results on the Traveling Salesman Polytope, an example of a complicated convex 
body studied in combinatorial optimization. 



1. Introduction 

Let V be a finite-dimensional real vector space. A set B C V is called convex if 
for every two points x,y G B the interval [x, y] = {ax + (1 — a)y : < a < l} also 
lies in B. A convex set B C V is a convex body if B is compact with a non-empty 
interior. 

Various notions of complexity for convex bodies have been discussed from a 
variety of points of view, see, for example, [Sz06] for a recent survey. In this 
paper, we suggest the following general approach to the computational complexity 
of convex bodies. 

Let B C V be a <i-dimensional convex body. We would like to design an efficient 
algorithm for the following membership problem 



(1.1) The membership problem for B 
Input: A point x £ V. 

Output: "Yes" if x E B and "No" iix£B. 
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The algorithm in Problem 1.1 is, of course, determined by the convex body B. It 
is important to note that we do not count the time and resources spent on creating 
such an algorithm towards the computational complexity of B. Once the algorithm 
is created, only the time needed to answer the question whether a given point x G V 
belongs to B counts as the complexity of B. 

At this point, it is not important to us what model of computation we use. For 
example, we can identify V = JBL d and assume that a point x = (£1,... , £d) is 
given by its real coordinates. The algorithm employs arithmetic operations and 
its complexity is the number of operations used (the real model). Alternatively, 
we may assume that the point x has rational coordinates and count bit operations 
instead (the bit model). We are interested in how the complexity of the algorithm 
depends on the dimension d. 

The algorithmic theory of convexity, see [G+93] , implies that as soon as there is a 
polynomial time algorithm in the Membership Problem 1.1 for a convex body B (in 
the bit model), there is a polynomial time algorithm in the problem of optimizing 
a given linear function on B (again, in the bit model). Hard optimization problems 
supply a variety of convex bodies for which the membership problem is likely to be 
hard. We concentrate on one model example, the (symmetric) Traveling Salesman 
Polytope. 

(1.2) Example: the Traveling Salesman Polytope. Let us fix an integer 
n > 4. In the space Mat n of nxn real matrices, we define TSP n as the convex hull of 
the adjacency matrices of the Hamiltonian cycles in a complete graph on n vertices. 
That is, for any Hamiltonian cycle C (a cycle that visits every vertex exactly once) 
of a complete undirected graph on the vertices {1, . . . , n}, we introduce the nxn 
matrix x = x(C) where 



and take the convex hull of all such matrices x(C). One can observe that TSP n 
has (n — l)!/2 vertices and that dimTSP n = (n 2 — 3n)/2, see Chapter 58 of [Sc03]. 
We define the ambient space V n as the affine hull of TSP n with the origin at the 
center a = ((%•), where 



The space V n consists of the nxn symmetric matrices with zeros on the diagonal 
and the row and column sums equal to 2. 

The facial structure of TSP n was studied in many papers but still remains a 
mystery, see Chapter 58 of [Sc03] for a survey. In some sense, unless the complexity 
hierarchy collapses, the facets of TSP n cannot be described. More precisely, there 
can be no algorithm of a polynomial in n complexity, which, given an inequality in 
V n decides whether it defines a facet of TSP n , see Chapter 4 of [Sc03]. 




1 if ij is an edge of C 
otherwise 
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Since the Membership Problem 1.1 can be quite hard, we also consider the 
computational complexity of the approximation problem: we want to find a set 
X C V approximating B such that the membership problem for X can be solved 
within a given complexity, preferably in time polynomial in dim V. 

(1.3) Measuring the quality of approximation 

In many important cases, the convex body B is symmetric, that is B = —B and 
we measure how well X approximates B by a number a > 1 such that 

X CB CaX. 

Often, B is not symmetric but has a natural center nevertheless, as in the case of 
the Traveling Salesman Polytope. In such a case, taking the center of B for the 
origin, we measure the quality of approximation in the same way. In many cases, 
but not in all, the set X is convex as well. 

This paper is meant to be a survey although it contains some new results, mostly 
in Section 5 and in Theorems 4.6 and 6.2. The paper is structured as follows. 

In Section 2, we discuss how well a symmetric convex body can be approximated 
by an ellipsoid, or, more generally, by an algebraic hypersurface of a given degree. 
In particular, for any e > and any symmetric convex body B we obtain a set X 
such that IcBc ctX for a = eVdimS and the membership problem for X can 
be solved in time polynomial in dim 5. 

In Section 3, we review results on approximation of convex bodies by polytopes 
with a controlled number of vertices or facets. We show that to approximate a 
(i-dimensional symmetric convex body within a constant factor, it is sufficient to 
take a polytope with vertices. The exponential in d number of vertices is also 
necessary in the worst shows the example of the Euclidean ball. 

In Section 4, we discuss approximations of convex bodies by projections of poly- 
topes with a controlled number of facets. Note that if X is a projection of a polytope 
with N facets then the membership problem for X is a linear programming prob- 
lem that can be solved in time polynomial in iV (in the bit model), cf. [G+93]. 
On the other hand, the number of facets of a polytope can grow exponentially 
if the operation of projection is applied, so potentially we get more flexibility for 
approximation. We describe several interesting phenomena here. First, we show 
that to obtain a good approximation of a symmetric convex body by a projection 
of a polytope with not too many facets, we may have to require the polytope to be 
non-symmetric. This is the case, for example, when the body is the cross-polytope 
(octahedron) . Then we discuss an amazing approximation, constructed by A. Ben- 
Tal and A. Nemirovski [BN01] of the Euclidean ball in IR d within a factor of (1 + e) 
by the projection of a polytope with only O (<ilne _1 ) facets. Finally, we discuss a 
general way to approximate an arbitrary convex body by a projection of a polytope 
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with a controlled number of facets. As an illustration, we state the result of the 
second author that for any e > the Traveling Salesman Polytope TSP n can be 
approximated within a factor of en by the projection of a polytope with n ^ 1 '*' 
facets. 

In Section 5, we present a construction which, for any <i-dimensional convex 
body B C V "almost approximates" its polar B° C V* by a section of a polytope 
with not more than e O(Vdind) ver ti C es. We consider V* embedded into the space 
C(B) of continuous functions on B, construct a polytope R C C (B) with e O(.Vd]nd) 
vertices and show that we obtain a good approximation of B° if we slightly "bend" 
V* C C(B) and intersect in with R. 

In Section 6, we discuss approximations by a section of the cone of positive 
semidefmite quadratic forms. Such approximations also lead to efficient algorithms 
for the membership problem. We review by now classic results on approximations 
of the cut polytope and its variations. We also present a general construction and 
illustrate it by the result of the second author that the polar TSP° of the Traveling 
Salesman Polytope can be approximated within a factor of en by a section of the 
cone of positive semidefmite quadratic forms in n ^ 1 ^^ variables. 

2. Approximation by algebraic hypersurfaces 

An ellipsoid E C V is a set defined as 

E={veV: q(v-v )<l}, 

where q : V — ► M. is a positive definite quadratic form and vq G V is a particular 
point called the center of the ellipsoid. As is known, see for example, [Ba97], for 
any convex body B C V there is a unique ellipsoid E C B that has the largest 
volume among all ellipsoids contained in B. It is called the John ellipsoid of B and 
it satisfies 

E C B C dE, 

where d = dimV and the origin of V is moved to the center of E. 

If B is symmetric then the John ellipsoid E is necessarily centered at the origin 
and satisfies 

E C B C VdE. 

For the Traveling Salesman Polytope TSP n the ellipsoid E is centered at the center 
of TSP n , touches the facets of TSP n defined by the equations Xij = and satisfies 

EGTSP n C ( n -*)^ E fo r n > 5 , 

see [BB05]. 

Since an ellipsoid E is defined by one quadratic inequality, the membership 
problem for E can be solved in O (<i 2 ) time for d = dimV . One can ask whether 
better bounds can be obtained by using polynomial inequalities of higher order. 

For symmetric convex bodies such a result was obtained in [Ba03]. We state 
some general prerequisites first. 
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(2.1) Symmetric convex bodies and norms. Polarity. With a symmetric 
convex body B C V one naturally associates a norm || • || B on V defined by 



\v\\ B = inf|A > : v G Asj. 



Approximating £? by an efficiently computable set X is equivalent to approximating 
|| • \\b by an efficiently computable function /. 

Let V* be the space of all linear functions £ : V — ► R. We recall that the polar 
B° C 7* of a convex body B C V is defined by 



= |^ G V* : ^(u) < f for all u G s}. 



The standard duality result states that (B°)° = B if B is a closed convex set 
containing the origin. 

The following result was proved in [Ba03] . 

(2.2) Theorem. For any symmetric convex body B C V and any integer k > 1 
there exists a homogeneous polynomial p : V — ► K. snc/i t/iat p is the sum of squares 
of homogeneous polynomials of degree k and 

P 1/2k (v) < \\v\\ B < a(d, k^^iv) for all v G B, 

where 

d + k-1^ 1/2k 



a(d,k) = y ^ j and d = dimB. 

In other words, the set 

X = {vEV: p(v) < l} 

approximates B within a factor of a(d, k). Note that the membership problem for 
X can be solved in O (d 2k ) time. If d 3> k 3> f , applying Stirling's formula, we get 

a(d, k) 




(k\) 1/2k 

It follows then that for any fixed e > we can choose k = k(e) large enough so that 
X approximates B within a factor of eV dim I? and the membership problem for X 
can be solved in time polynomial in dim 5. It is not clear whether the set X can 
always be chosen convex, although this is the case in many situations. 
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Sketch of proof of Theorem 2.2. Let V* be the dual space of all linear functionals 
£ : V — >■ R and let 



B° = {I e V* : t(v) < 1 for all v E B} 
be the polar of B. By the standard duality argument, we can write 

\\v\\ B = max|£(v) : ^B°). 
Let V m be the tensor product of k copies of V, so (V® k )* = (V*f k . Then 

= max<J> fc : IgB°}. 

Now, let us define a symmetric convex set C C (y*^® k as the convex hull 

C = conv(l m ,-l m : £eB°y 

Then 



Ml r = max 



Finally, we approximate C by an ellipsoid E C C The crucial consideration is that 
since C lies in the symmetric part of (y*)® k , we have 



dimC < 



/dimV + k-1 



and hence E approximates C within a factor of a k (d, k). 
It remains to notice that the formula 



V^y = max|L(t;® fc ) : L G i?| 



indeed defines a polynomial p : V — >■ K. of degree 2/c which is a sum of squares of 
polynomials of degree k. □ 

There are several open questions related to Theorem 2.2. 

(2.3) The best approximation? 

Suppose that we want to approximate an arbitrary norm \\v\\b by an expression 
p 1 / 2fc (v), where p : V — ► f is a homogeneous polynomial of degree 2k. Is the 
coefficient a(d, k) of Theorem 2.2 best that we can get? This is so in the case of 
k = 1, since the <i-dimensional cube 

Jd={(£i,- -.,&): 161 <1 for i = l,...,d} 



or the (i-dimensional cross-polytope (octahedron) 



d 

Od={(£i,...,&): El^l^ 1 } 

i=i 

cannot be approximated by an ellipsoid better than within a factor of y/d, cf . [Ba97] . 
Similarly, in the non-symmetric situation, the (i-dimensional simplex 

d+l 

A d = {(£1,... : = 1 and &>0 for i = l,...,d + l} 

i=l 

cannot be approximated by an ellipsoid better than within a factor of d. 

For k > 1 it is not clear whether the bound of Theorem 2.2 is optimal. One can 
show, however, that the octahedron Od cannot be approximated by a hypersurface 
of degree 4 or 6 better than within a certain factor cy/d for some absolute constant 
c > [Ba03]. The idea of the proof is as follows. Once a convex body B and the 
degree 2k are fixed, there exists the best approximating polynomial p. Furthermore, 
the polynomial p can be chosen to be invariant under the group of the symmetries of 
B. In the case of O d , it follows that one can choose p to be a symmetric polynomial 
in £i, . . . , from which the estimate can be deduced. 

It is interesting to note that while the cube Id and the octahedron Od demon- 
strate similar behavior with respect to ellipsoidal approximations, their behavior 
with respect to higher degree approximations very much differ. Indeed, the cube Id 
can be approximated by a hypersurface of degree 2k within a factor of d 1 ^ 2h . For 
that, one can choose 

d 

i=i 

(2.4) Convex approximation? 

Although in many special cases the function p 1 / 2k (-) constructed in Theorem 2.2 
turns out to be a norm, there is no reason to believe that it is always a norm (but 
there are no explicit counterexamples either). Equivalently, there is no apparent 
reason why the set 

X = {veV: p(v)<l} 

should be convex. It is an interesting question whether one can choose p to be a 
norm and still get the same order of approximation. One natural candidate would 
be 

(2.4.1) p(v) = [ £ 2k (v) dixit), 

J B° 
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where \x is a certain probability measure on the polar B° . One can show that for 
the normalized Lebesgue measure \i the bounds generally are much weaker than 
those of Theorem 2.2. A natural candidate is the exterior angle measure defined as 
follows. 

Let S d ~ 1 C R d be the unit sphere in Euclidean space R d endowed with the scalar 
product () and let v be the Haar probability measure on Let K C M d be a 

convex body. For a closed subset X C K, let us define 

uk(X) — v\c£l S d_1 : max(c, x) = max(c, x) >. 

I xex xek J 

In words: the measure of a closed subset X is the proportion of linear functions of 
length 1 that attain their maximum on K at a point x G X. Clearly, the probability 
measure \ik depends on the convex body K as well as on the Euclidean structure 
in R d . 

Although no bounds have been proved, it seems plausible that by choosing a 
scalar product in V in which the John ellipsoid E of B is the unit ball and letting /j, = 
Hb° in (2.4.1), one obtains a norm p 1 /^ which gives a similar order approximation 
as in Theorem 2.1. 

3. Approximations by polytopes 

Any convex body B C V can be arbitrarily well approximated by the convex hull 
of a sufficiently large finite subset X C B. Consequently, the membership problem 
for B is replaced by the membership problem for the polytope P = conv(A) whose 
complexity (in the bit model) is polynomial in the cardinality \X\ of A. In the dual 
setting, any convex body can be arbitrarily well approximated by the intersection 
of a sufficiently large set X of halfspaces. Then the complexity of the membership 
problem is linear in \X\ in both the real and the bit models. How large should 
X C B be so that conv(A) approximates B reasonably well? 

The following result is well-known, see for example, Lemma 4.10 of [Pi89]. 

(3.1) Lemma. Let B C V be a symmetric convex body. Then, for any 1 > e > 

there is a subset X C B such that 

\X\ < ( 1+< ^) f° r d = dimV 

and 

PCBC- -P for P = conv(A). 

Proof. Let || • \\b be the norm associated with B, see Section 2.1. Let X C B be a 
maximal (under inclusion) subset such that 



xi — X2\\b > e for all distinct x±,X2 G X. 
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Since X is maximal, for every x G B there is a point y G X and a point z E eB 
such that x — y = z. In other words, 

BcI + eBcP + e5. 

Iterating, we get 

B cP + eP + e 2 P + ... + e k P + e k+1 B 
and taking the limit we conclude that 

B C — - — P. 

1 — e 

Next, we estimate the number of points in X. We notice that 

(x 1 + n (x 2 + = for distinct Xi,x% G X 



Furthermore, 



Hence 



from which 



U(* + ^)c(l + i)*. 



x6X 



vol (i + |)b> |X|vol(|s 



\x\ < 



2 + e^° 



as required. □ 

Although in many cases the bound can be slightly sharpened, the bottom line 
is that the bound on the number \X\ of points is exponential in the dimension. In 
fact, simple volume estimates show that even if B C M d is the Euclidean ball, the 
number of points needed to approximate B within a constant factor is exponentially 
large in d. 

The following result is proved, for example, in [Ba97]. 

(3.2) Theorem. Let B C M, d be the unit ball and let X C B be a set such that for 
P = conv(X) we have 

P C B C aP. 

Then 



|X|£eXp \2^ 

Similarly, if P is the intersection of a set X of half spaces such that 

P C B C aP 

then 

\X\ > exp 
9 
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4. Approximation by projections 



Let B C V be a convex body. Suppose that we manage to construct a poly tope 
P C W, where W is some other vector space, and a linear transformation T : 
W — > V such that the image Q = T(P) approximates B reasonably well. Given an 
x EV , testing whether igQ reduces to testing whether the affine subspace T~ 1 (x) 
has a non-empty intersection with P, which is a linear programming problem. In 
particular, if P is defined by N linear inequalities, the complexity of the membership 
problem for Q is bounded by a polynomial in N (in the bit model). On the other 
hand, the number of facets of Q can be exponentially large in N, which shows 
that, in principle, Q may provide a fairly good approximation even to rather "non- 
polytopal" convex bodies B. 

In the dual setting, we are interested in approximating the body B by a section 
of a polytope. Namely, we want to construct a vector space W D V and a polytope 
P C W such that the intersection Q = P n V approximates B reasonably well. 
If P is defined as the convex hull of N vertices, the membership problem for P, 
being a linear programming problem, has the complexity that is polynomial in N 
(in the bit model). On the other hand, the number of vertices of Q = P n V can 
be exponentially large in N, which, in principle, can lead to good approximations 
of B. 

One can observe that a polytope Q is a projection of a polytope with at most 
N facets if and only if Q is a section of a polytope with at most N vertices, so the 
two approaches are essentially identical. 

(4.1) Good sections are good projections and vice versa. 

Let V C W be a pair of spaces, let P C W be a polytope with N vertices 
and let Q = P n V. Since P has N vertices, there exists a simplex A C M. N , see 
Section 2.3, and a linear transformation T : M N — > W such that P = T(A). Then 
U = T _1 (l/) C 1^ is a subspace, the polytope P' = A n U has at most N facets 
and Q = T(P'). In other words, a polytope that is a section of a polytope with at 
most iV vertices can be represented as a projection of a polytope with at most iV 
facets. 

Vice versa, let P C W be a polytope with iV facets, let T : W — > V be a 
linear transformation and let Q = T(P). Fixing a scalar product in W, we identify 
V with a subspace V C W and T with the orthogonal projection of W onto V. 
Assuming that P contains the origin in its interior, let us consider the polars P° 
and Q°. One can see that Q° = P° fl V. Since P° is a polytope with iV vertices, 
by the above reasoning Q° can be represented as a projection of a polytope with 
at most iV facets. Dualizing again and using that (<2°)° = Q, we conclude that Q 
can be represented as a section of a polytope with at most iV vertices. 

Approximations by projections are well suited for taking intersections and direct 
products. 

(4.2) Operations on projections. 
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For i = 1,2, let Pi C Wi be polyhedra, let T; : Wi — > V be linear transforma- 
tions, and let Qi = Tj (Pi). Let Q = (Q\ fl Q 2 ) C V. We observe that Q is the 
image of the polyhedron P C Wi © W^, 

P=|(xi,x 2 ): xi G Pi, x 2 G P2, and Ti(xi) = T 2 (x 2 )| 

under the linear transformation defined by T(x\,X2) = Ti(x\). If Pj has at most 
iVj facets for z = 1, 2 then P has at most N = N\ + N 2 facets. 

For i = 1, 2, let Pj C Wi be polyhedra, let Tj : Wi — > V$ be linear transforma- 
tions, and let Qi = Ti (Pi). Let Q = (Q\ x Q 2 ) C V\ © V2. Then Q is the image of 
the polyhedron P = (Pi x P 2 ) C W\ © under the linear transformation defined 
by T (x\, X2) = (T\ (x\) ,T 2 (x 2 )). Similarly, if Pj has at most Ni facets for i = 1, 2 
then P has at most N = Ni + N 2 facets. 

One interesting feature of approximations by sections and projections is that 
it breaks symmetry. Namely, to obtain a sufficiently close approximation of a 
symmetric convex body B by a projection of a polytope P with not too many 
facets, we may have to choose the polytope P to be non- symmetric. 

Let Oa C M d be the standard octahedron (cross-polytope) , see Section 2.3. 
Clearly, Oa has 2d vertices and hence there is a simplex A 2 ^_i C M. 2d with 2d 
vertices and a linear transformation T : R 2d — > M. d such that T (A 2 d-i) = Oa- In 
other words, can be represented as the projection of a polytope with 2d facets. 

Suppose, however, that we want to construct a symmetric polytope P C W 
and a linear transformation T : W — ► M d such that Q = T(P) approximates Oa 
within a factor of 2. One can show that the number iV of facets of P has to be 
exponentially large in d. 

(4.3) Theorem. Let Oa C M. d be the cross-polytope 

a 

i=i 

and suppose that P C W is a symmetric polytope with N facets and T : W — ► M d 
is a linear transformation such that 

QcO d c2Q for Q = T(P). 

Then 

N > e cd for some absolute constant c > 0. 

Proof. The proof uses the notion of the type 2 constant of a Banach space, see 
[Pi89] and [To89]. 

Let V be a finite-dimensional vector space, let B C V be a symmetric convex 
body and let || • ||_b be the corresponding norm. The type 2 constant of B is the 
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smallest number k = k(B) > such that for any set of vectors xi, . . . , x m E V we 
have 

2 



E J^eiXi < K 2 (B)J2\\ x i\\ 2 B, 



i=i i=i 
where the expectation is taken with respect to independent random signs e^: 

1 with probability \ 
-1 with probability |. 

One can observe that k(B) = 1 if B is the unit ball in some Euclidean metric and 
that if T : V — ► W is an invertible linear transformation then k(B) = k (T(B)) for 
any symmetric convex body B G V. Since every convex body can be approximated 
by an ellipsoid, it follows that k(B) exists (in particular, is finite) for all symmetric 
convex bodies B. 

Let B C W be a symmetric convex body and let U C W be a subspace. It is 
immediate that 



(4.3.1) 



k(BHU)< k(B). 



Furthermore, suppose that T : W — > V is a linear surjection. Then 
(4.3.2) k(T(B)) < k(B). 

To prove (4.3.2), we notice that 

\\Ty\\ T {B) < \\y\\ B for all y E W 

and that 

for every x E V there is y E W such that T(y) = x and \\y\\B = ||^||t(b)- 

To establish (4.3.2), let us choose any x\, . . . , x m E V and then yi, . . . , y m E W 
such that 

T(yi) = Xi and ||yj||s = ||^||t(b) for i = l,...,m. 



Then 

m 

E IIE 
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i=l 



T(B) 



<v\\j2e iyi < K 2 (B)J2\\yi\\l = K 2 (B)J2 



\T(B) 



i=i 



i=l 



i=l 



and the proof of (4.3.2) follows. 
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Suppose now that P is a symmetric poly tope with N facets. Then, up to a linear 
transformation, P can be represented as a section of an iV-dimensional cube In, 
see Section 2.3. Hence k(P) < k (In) and if Q is a projection of P, we have 

k(Q) < k(P) <k(I n ). 

Furthermore, under the conditions of the theorem, we have 

K(O d ) < 2k(Q) < 2k (I n ). 

On the other hand, one can show that 

K(O d ) > Vd and k(I n ) < cVlniV 

for some absolute constant c > 0, see Section 4 of [To89]. This proves that N has 
to be exponentially large in d. □ 

Theorem 3.2 shows that the <i-dimensional Euclidean ball is not approximated 
very well by a polytope with a subexponential in d number of vertices or facets. 
The situation changes radically if we consider approximations by a projection of a 
polytope with a small number of facets or by a section of a polytope with a small 
number of vertices. For example, the intersection of a 2<i-dimensional octahedron 
Oid with a random <i-dimensional subspace of M? d approximates the <i-dimensional 
Euclidean ball within an absolute constant, see, for example, [Ba97]. The results of 
[B+89] on approximation of zonoids by zonotopes imply that for any d and e > 
there is an N = O* (de~ 2 ) , where * stands for some logarithmic factors, and a linear 
transformation T : R N — > R d , such that the image T (I N ) of the cube I N C R N 
approximates the unit ball B C M. d within a factor of (1 + e). Note that in this 
case we approximate a symmetric convex body by the projection of a symmetric 
polytope. 

A. Ben-Tal and A. Nemirovski [BN01] provide a remarkable construction of a 
polytope P with N = O (dine -1 ) of facets whose projection T(P) approximates 
Euclidean ball B C W 1 within a factor of (1 + e). We sketch the construction below. 
Unlike in the case of zonotopal approximation of [B+89], the polytope P is not 
symmetric. 

(4.4) A tight approximation of the Euclidean ball by a projection. Let 

d 
i=l 

be the <i-dimensional unit ball. The first idea of [BN01] is to reduce the general 
case to that of d = 2. Let us consider the (d + 1) -dimensional round cone 

d 

C d ={(£i,...,£ d ;r): ^4 2 <r 2 and r > o}. 

i=i 
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One can observe that a close approximation of by a projection of a polyhedron 
with at most N facets results in a close approximation of by a projection of a 
polytope with at most N facets and vice versa. 

Let us write d = r + b for some positive integers r and b and let us consider the 
corresponding round cones of lower dimensions: 



CV = {(£i,...,£r;p): J2^<P 2 and P> } and 

r 

<7 6 = {(t/i,..., 7/ 6 ;0): $> 2 </3 2 and /3 > o}. 



i=i 



In the space M d+3 with the coordinates £i, . . . , £5; ?7i> • • • , ?7r; P, /3, t we consider the 
subset 

x = I (£1,... ,£ r ;vi,---Vb;p,P,r) : 

$^ 2 <P 2 , J>?</? 2 , p 2 + /3 2 <r 2 and 

i=l i=l 
P,P,T >()}. 

Clearly, Ca is the image of X under the projection which forgets p and (3. On the 
other hand, we have 

X = (C r X C b X Vr) n X V^J , 

where 

C2 = {(p,/3,r) : /3 2 + p 2 <r 2 and p,/3,r>o} 

is the "quarter" of the 3-dimensional round cone and V T and V^ v are the appropriate 
coordinate subspaces. 

Using the operations of direct product and intersection, see Section 4.2, one can 
show that good approximations of C r , C\,, and C2 by projections of polyhedra leads 
to a good approximation of Ca by a projection of a polyhedron. 

Essentially, the problem boils down to an efficient approximation of the quarter 
of the disc: 

^={(£,77): Z 2 + V 2 <1 and £,r/>0}. 

For that, A. Ben-Tal and A. Nemirovski provide the following ingenious construc- 
tion. 

Let us define the sequence of transformations R n of R 2 by 

f = f cos h 77 sin — 



. 7T 7T 

— t sin h 77 cos — 
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Geometrically, R n is a clockwise rotation through an angle of n/2 n followed by 
the reflection in the £-axis if the obtained point lies in the lower halfplane. If we 
pick a point z = (£,77) with £,77 > and apply the sequence of transformations 
i?2, -R3, -R4, • • • , then the resulting sequence of points has its limit on the interval 
< £ < 1, 77 = if and only if z G ~B~ 2 . 

Let us choose a positive integer m. In the space of variables ik,i]k for k = 
1, . . . , m we define the polyhedron P m by the equations and inequalities 

£ fc = £ fc _i cos ^ + ?7 fc _i sin 

77fc > sin + 77 fc -i cos ^ 

77fc > Cfc-i sin J - ?7 fc -i cos J for k = 2,... ,m 



and 



7T 

Ci,77i>0, 0<£ m <l, r/ m < — . 



The projection of P m onto the coordinates (£1, 771) approximates B 2 within an error 
exponentially small in m. 

(4.5) A general approximation construction. Let B C M d be a convex body 
we want to approximate by a projection of a polytope with at most N facets, or, 
equivalently, by a section of a polytope with at most iV vertices. Here is a general 
construction. 

Without loss of generality, we assume that B contains the origin in its interior and 
hence we may view B as the polar of its own polar B° C M, d . Next, we approximate 
B° by a sufficiently dense finite subset X C 5° so that 1° approximates £? well 
enough. 

Let us consider the space R x of all functions / : X — > R and let A C R x be 
the affine subspace of all affine functions / : X — > R whose average value on X 
is 1. In other words, A consists of the functions of the type f(x) = (c,x) + a for 
some c G R d and a G R such that, additionally, 

m ^ ,{x) = L 

1 1 xex 

It is not hard to see that X° can be viewed as 

X° = |/ G A : f(x) > for all x G X j. 

For a subset Y - C X let 5y G R x be the indicator of Y, that is 



S Y (x) 



1 if x G Y 

if x ^ Y. 
15 



We write simply 5 X instead of 5{ x y. Now we can write 

X° = A n co (S x : x G X) , 

where "co" stands for the conic hull of the set, that is, the set of all non-negative 
linear combinations of the elements of the set. 

Let us choose a family T C 2 X . We define the approximation of Bp as 

B T = Anco(5 F : F e J 7 ) . 

Clearly, Bp is the intersection of an affine subspace with a polyhedron with at most 
\J-\ vertices. One can expect that the finer T gets, the better approximation Bjr 
to B we obtain. For example, we may choose X as the intersection of B with a 

sufficiently dense grid, X = (b n (eZ) d ) and then choose consisting of the sets 

of points lying on an affine coordinate subspace of codimension k. Generally, we 
will have \F k \ = O (d k e~ k ). 

Using this approach, the second author was able to obtain the following approx- 
imation result for the Traveling Salesman Polytope TSP n , see Example 1.2. 

(4.6) Theorem. Let us choose an e > and let V n be the ambient space ofTSP n . 
Then there exists a polytope P n C W n with at most N = O (n 4 / e ) facets and a 
linear transformation T : W n — ► V n such that for Q n = T(P n ) one has 

Q n C TSP n C (en)Q n . 

The proof of Theorem 4.6 is rather technical and will be presented elsewhere. 

(4.7) Remarks and open questions. 

Given a series of combinatorially defined polytopes Q n , such as the Traveling 
Salesman Polytope TSP n , to construct a simpler polytope P n whose projection 
approximates Q n , by now is a well-established technique of "lift and project" in 
combinatorial optimization, see, for example, [BaOl]. In particular, a remarkable 
success was achieved in showing that certain exponentially large (in ri) families of 
facets of certain series of polytopes Q n can be obtained as projections of only poly- 
nomial in n families of facets of P n , see [LS91]. However, not much is known about 
how well the projections approximate metrically, even though various techniques 
were compared combinatorially, see [La03]. 

We also note that the authors of [BN01], while constructing their remarkable 
approximation of the Euclidean ball by the projection of a polytope with a small 
number of facets (cf. Section 4.4), were motivated by very practical questions. 
Namely, they used their approximation to reduce convex quadratic programming 
problems to linear programming problems (see [BN01] for details) and hence use 
the linear programming solver available to them to solve quadratic programs. 

Some basic questions regarding approximations by projections remain unan- 
swered. 
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(4.7.1) Obstructions to being a projection. Let Q C V be a polytope. How can 
one possibly prove that Q cannot be the projection of a polytope P with at most 
iV facets? In the case when Q is symmetric and P is required to be symmetric, a 
possible argument goes via the type 2 constant, see Theorem 4.3. Another example 
is the proof of [Ya91] that TSP n cannot be a projection of a polytope P n with a 
polynomial in n number of facets provided the projection respects the symmetries 
of TSP n . In general, if the lifting P is not required to be symmetric, no viable 
argument seems to be known. It would be interesting to find out whether some 
appropriate notion of a non-symmetric type can be of help. 

(4.7.2) The quality of approximation of a general body. How well can a general 
(i-dimensional convex body B be approximated by the projection of a polytope 
with at most N facets? In particular, the following question seems to be of interest. 
Suppose that B is symmetric and that the projection of a polytope P approximates 
B within a factor of 2. Is it true that in the worst case the number N of facets 
of P should be at least exponential in d, N > c d for some absolute constant c > 
1? In Section 5 we discuss a certain construction which suggests that maybe the 

exponential bound can be broken and that we can have N = O* ^c^j , where, as 

usual, * stands for some logarithmic factors. 

(4.7.3) Approximation of the V ball. Let p > 1 and let 

d 

5(d,p) ={&,...,&): ££f<l} 

i=i 

be the unit ball in the V norm. One can observe that the first step of the construc- 
tion of Section 4.4 can be extended to B(d,p) thus providing an approximation 
of B(d,p) within a factor of (1 + e) by the projection of a polytope with at most 
N = O (rfe -1 ) facets. It is not clear whether for a general p one can replace e _1 by 
lne -1 , though this is definitely the case for p = 1, 2, and +oo. 

5. A "soft" approximation of a symmetric convex body 

Let B C V be a symmetric convex body. We identify V* with the subspace of 
linear functions in the space C (B) of all continuous functions / : B — ► HL Then 
B° can be identified with the set of linear functions / : B — > R such that f{x) < 1 
for all x £ B. 

In this section, we prove the following main result. 

(5.1) Theorem. Let B C V be a symmetric d-dimensional convex body and let 
C (B) be the space of all continuous functions on B. Then there exists a polytope 
R C C (B) such that the following holds. 

(1) For all h G R we have 



h(x) < 1 for all x G B; 
17 



(2) For any £ G B° there exists a function h G R such that 
\£(x) - h(x) | < <y£ 2 (x) for all xeB 



and some absolute constant 7. 
(3) The polytope R has at most exp jav^i In vertices for some absolute con- 
stant a. 

If we could claim that the function h in Part (2) of the theorem is linear, then 
we must have had h(x) — £(x) in Part (2) and we would have obtained the repre- 
sentation 

B° = R n V* 

of B° as the intersection of a polytope with at most exp |av^ln<i| vertices and 

a subspace V* . By duality, that would have implied that B is the projection of a 

polytope with at most exp jav^i hi <ij facets. In general, however, h is not a linear 

function, but, as will follow from the proof, is a piecewise polynomial. If £ G eB° 
for some < e < 1 then h approximates £ within an error of O (e 2 ) , so the points 
of B° that are closer to the origin are better approximated. Intuitively, we obtain 
a set close to B° if we slightly "bend" V* and then intersect it with R. 

Proof of Theorem 5.1. Since B is symmetric, we can find an ellipsoid E D B 
centered at the origin such that 

-^=E C B C E. 

Vd 

On the other hand, approximating the ellipsoid by the projection of a polytope (see 
Section 4.4), one can construct a polytope P C W with iV = 0(d) facets and a 
linear transformation T : W — >■ V such that 

-T(P) C E C T(P). 
2 

Summarizing, 

(5.1.1) ^ff {P) CBC T{P) - 

Without loss of generality, we assume that P C W is full-dimensional and contains 
the origin in its interior. Suppose that 

P={weW: gi(w) < 1, i = 1,... ,ivj, where g l : W — >R 
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are linear functions. Let us choose the smallest positive integer k > 2y/d. For a 
multiset I of numbers from 1 to N of cardinality at most k (counting multiplicities), 
we let 

gi = i-Y[(i- gi ). 

iei 

Hence gi : W — > R are polynomials, deg gi < k. It is immediate that gi(w) < 1 
for all w E P. 

Suppose now that I E B° , so I : B — > R is a linear function such that £(x) < 1 
for all x E B. In view of (5.1.1), we have i(x) < k for all x G T(P). Let / : W — > R 
be the lifting of I defined by f(w) = I (T(w)) for w G W. Hence / : W — > R is a 
linear function and f(w) < k for all w G P. Therefore, k~ 1 f(w) < 1 for all w G P 
and hence 

AT 1 / G conv(0, gi : z = l,...,iV). 

Therefore, 

(5.1.2) F = 1 - (1 - k- x f) k G conv^O, gi : |/|<ife). 

Now, since B symmetric, we have \£(x)\ < 1 for all x G B. Therefore, for all w G P 
such that T(w) G B, we have \f(w) \ < 1 and hence 

(5.1.3) \F(w) - f(w)\ < -ff 2 (w) for all weP such that T(w) G B. 

Now we are ready to define R C C(B). Let us fix a scalar product in W and 
hence the Lebesgue measure on every affine subspace of W. Let us define hi G C(B) 

hi(x) = the average value of gi(w) over all w E P (~)T 1 (x). 

Let us define R to be the convex hull of the origin and all the functions hi as 
/ ranges over all multisets with the elements from {1, . . . ,N} and of cardinality 
at most /c, counting multiplicities. Clearly, the number of vertices of R does not 
exceed N h , so Part (3) follows. 

Since hi(x) < 1 for all / and all x E B, Part (1) follows as well. 

To prove Part (2) we choose 

h(x) = the average value of F(w) over all w E P H T _1 (x), 
where F is defined by (5.1.2). Clearly, h E R and Part (2) follows by (5.1.3). □ 

Let I E V* be a linear function I : B — > R. There seems to be no efficient way 
to check whether there is a function h E R such that \£(x) — h(x)\ < 7^ 2 (x) for 
all a; G B. We can relax the condition by replacing the uniform distance by the 
distance in the L 2 (B, fi) norm for some Borel probability measure \i on B: 

11/1,2 = (X /2 ^ 
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Then, checking whether for a given <6 7* there is an h G R such that 

\\£-h\\ 2 <M\Z 2 h 

becomes a problem of convex quadratic programming which can be solved roughly 

in exp |o(v / rfln<i) j time (in the bit model). 

Let us choose a sufficiently small e > 0. Let us "accept" a given linear function 
£ G V* if there exists an h G R such that 

P-/*|| 2 <7«2 

and "reject" it otherwise. Hence we get an algorithm of exp com- 
plexity such that given a linear function £ : B — > R 

(i) the algorithm accepts £ if £(x) < e for all x G -B; 

(ii) if the algorithm accepts £ then there is a function h : B — > H. with < 1 
for all x e B and ||£ — h\\ 2 < 1 1 1 1 2 5 where 7 is an absolute constant. 

This is somewhat similar to the situation of the "property testing" in computer 
science, see, for example, [Go99]. 

One can observe that the estimates of Theorem 5.1 can be extended in a number 
of ways. If we know that B can be approximated by an ellipsoid within a factor 
p > 1 then we can construct a polytope R with the number of vertices not exceeding 
exp {otp In d}. More generally, we can require that R has not more than exp {ak In d} 
vertices for k < p if we replace the condition (2) in Theorem 5.1 by \kp~ 1 £(x) — 
h(x)\<-fk 2 p- 2 £ 2 (x). 

6. Approximations by a section of the cone 
of positive semidefinite quadratic forms 

Let W be a finite-dimensional real vector space and let S(W) be the space of all 
quadratic forms q : W — >■ R. Let S + (W) C S(W) be the convex cone of positive 
semidefinite quadratic forms, that is, the quadratic forms such that q(w) > for 
all w G W. The membership problem for the cone S + (W) can be solved in time 
polynomial in dimVF both in the algebraic and the bit models of computation. 
Roughly, checking that q is positive semidefinite can be done by a certain variation 
of the Sylvester criterion. Identifying W = M d , we conclude that q is positive 
semidefinite if and only if the form 

qe(x)=q(x)+e(e i + ... + e d ) 

is strictly positive definite for all e > 0. Sylvester's criterion then implies that the 
the d principle minors of q e should be positive for all e > 0, which, in turn, reduces 
to checking that d univariate polynomials in e are positive for all e > 0. 

Let B C V be a ^-dimensional convex body. We may try to approximate B by a 
set X which is the intersection of the cone Ss r (W) for some W and a <i-dimensional 
affine subspace identified with V. The main accomplishments of this approach are 
associated with approximations of the cut polytope, see [DL97] . 
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(6.1) Approximating the cut polytope. For n-vectors x = (£1,... , £ n ) and 
V = • • • ? ^n), let us define the nxn matrix as the matrix with the (i, j)th 
entry equal to ^rjj and let Mat n be the vector space ofnxn matrices. 

Let us identify the space S (M n ) of quadratic forms on MJ 1 with the space ofnxn 
symmetric matrices. The cut polytope CUT n C S (R n ) is defined as the convex hull 
of all n x n matrices x <g> x for all vectors x = (£i, . . . , £ n ) with £j = ±1. As is the 
case with the Traveling Salesman Polytope, the membership problem for the cut 
polytope is NP-complete, cf. [DL97]. 

We consider CUT n as a subset of the space 5" (M n ) of symmetric nxn matrices. 
Let A C S (W 1 ) be the affine subspace consisting of the matrices with l's on the 
diagonal. It turns out that the intersection S + (M. n ) C\A approximates CUT n within 
a logarithmic factor with respect to the center at the identity matrix I £ S (M n ): 

CUT n C inS+(r) C c\nn(CUT n ) 

for some absolute constant c > 0. The logarithmic factor cannot be improved, see 
[A+06]. 

A variation of CUT n is what we call the asymmetric cut polytope ACUT n defined 
as the convex hull of all matrices x <8> y where x and y are n- vectors with the coordi- 
nates £i,77j = ±1. Again, the membership problem for ACUT n is NP-complete. It 
turns out that ACUT n can be tightly approximated by the projection of a section 
of a cone of positive semidefinite quadratic forms. Namely, let A C S (M 2n ) be 
the affine subspace of symmetric 2n x 2n matrices with l's on the diagonal. Let 
: S (R 2n ) — ► Mat n be the projection 

B x \-x 
c)- x 

and let 

Q n = <p(Ans(R 2n )). 

In words: Q n is the set of all possible nxn upper right corner submatrices of a 
2n x 2n positive semidefinite matrix with l's on the diagonal. It is not hard to see 
that ACUT n C Q n . Indeed, the matrix z ® z, where z = (£i, . . . , £ n , 771, . . . , r] n ) 
with i]j = ±1 is positive semidefinite with l's on the diagonal and has the matrix 
x ®y for x = (£1, . . . , £ n ) and y = (771, . . . , r] n ) as its upper right corner submatrix. 
It turns out that Q n approximates ACUT n within a constant factor: 

Q n C ACUT n C kQ u 

for some absolute constant k, called the Grothendieck constant. Its exact value is 
not known, but it is known that 

1.5708 w - < k < w 1.7822, 

2 ~ " 21n(l + v / 2) 
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see [A+06] and [AN06] for survey and recent developments. 

We note that the membership problem for the projection of a section of a cone 
of positive semidefinite quadratic forms is an instance of semidefinite programming, 
which can be solved in polynomial time (in the bit model), though only approxi- 
mately, cf. [K102]. 

Now we describe a construction for approximating a general convex body by a 
section of the cone of positive semidefinite forms. 

Let B C V be a convex body containing the origin in its interior. As in Section 
2, we think of B as the polar B = (B°)° to its own polar B° C V*. Hence 

B = {xeV: £{x)<l for all £ G 

Let us choose a positive integer k and let W k be the space of all polynomials 
p : V* — > M, degp < k. In particular, for any fixed k, the dimension of W k is 
bounded by a polynomial in the dimension of V. 

Let us choose a Borel probability measure \i on B° . For a point v G V, let 
Qv '■ Wk — > K. be the quadratic form defined by 

q v (p) = [ (1 - £(v))p 2 (£) dfi(£) for p G W k . 

J B° 

Clearly, if v G B then the form q v is positive definite and the set 

{q v : vev}cS(W k ) 

is an affine subspace. This allows us to define an approximation X k to B by 

X k = [veV: q v eS+(W k )}. 

Thus B C X k for all k. To show that X/» approximates B reasonably well, we have 
to show that for any point v ^ B which is sufficiently far away from B we can find a 
polynomial p : V* — > K. which takes large values on I G B° such that £(v) > 1 and 
small values everywhere else on B°. Besides, the value of fi{£ G B° : £(v) > l} 
should be sufficiently large, in particular, the exterior angle measure \ib° discussed 
in Section 2.4 can be of help. 

In the case of the Traveling Salesman Polytope, the second author obtained the 
following result. 

(6.2) Theorem. Let us choose an e > and let V n be the ambient space ofTSP n . 
Let TSP° C V* be the polar of the Traveling Salesman Polytope with respect to its 
center as the origin. Then, there exists a set Q n C V* isometric to the section of 
the cone of positive semidefinite quadratic forms in n 01 - 1 /^ variables by an affine 
subspace such that 

Q n C TSP: C (en)Q n . 

The proof is presented in [Ve06]. It is not clear whether the bound is sharp. 
For example, the vertices of TSP° corresponding to the standard facets Xij = of 
the TSP actually belong to Q n while the vertices of TSP° corresponding to the 
subtour elimination facets, see Chapter 58 of [Sc03], lie in a n Q n for a n = O (\/n). 
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